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Abstract. We consider random extended surface perturbations in the transverse 
field Ising model decaying as a power of the distance from the surface towards a pure 
bulk system. The decay may be linked either to the evolution of the couplings or 
to their probabilities. Using scaling arguments, we develop a relevance-irrelevance 
criterion for such perturbations. We study the probability distribution of the surface 
magnetization, its average and typical critical behaviour for marginal and relevant 
perturbations. According to analytical results, the surface magnetization follows 
a log-normal distribution and both the average and typical critical behaviours are 
characterized by power-law singularities with continuously varying exponents in the 
marginal case and essential singularities in the relevant case. For enhanced average 
local couplings, the transition becomes first order with a nonvanishing critical surface 
magnetization. This occurs above a positive threshold value of the perturbation 
amplitude in the marginal case. 



1. Introduction 

Quenched, i.e., time-independent disorder has a strong influence on the nature of 
quantum phase transitions which take place at zero temperature [1] . Many interesting 
features of the effect of randomness can be observed in one-dimensional (Id) systems 
for which several exact results have been obtained. Recently, the Id random 
transverse-field Ising model (TIM) has been the subject of much interest. It is defined 
by the Hamiltonian: 

ft = -^E( J ^+i + Vf), (i-i) 

I 

where af and af are Pauli matrices at site I and the exchange couplings, J l , the 
transverse fields, are random variables. 

For homogeneous independently distributed couplings and fields, Fisher [2, 3] 
obtained many new results about the static properties of the random TIM in (1.1) 
using a real space renormalization group method which is believed to become exact 
at large scales, i.e., sufficiently close to the critical point. Later, Fisher's conjectures 
have been checked numerically [4-6] and new results have been obtained, especially 
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about the dynamical properties of the model at the critical point [7, 8] as well as in 
the Griffiths phase [9,10]. 

In many physical situation the disorder is not homogeneous. A free surface or a 
defect in the bulk can evidently induce a disorder which is inhomogeneously distributed 
in its vicinity. In the present paper we shall consider a random version of the Hilhorst- 
van Leeuwen (HvL) model [11-13]. In this model, the couplings are modified near a 
free surface and the amplitude of the perturbation decays according to a power law 
l~ K where / denotes the distance from the surface. A relevance- irrelevance criterion 
has been proposed for such smoothly inhomogeneous perturbations, showing that 
marginal behaviour is obtained when the decay exponent k = \jv where v is the 
correlation length exponent of the pure system [14-16]. The critical behaviour is 
then characterized by local exponents which vary continuously with the perturbation 
amplitude. For a slower decay the perturbation becomes relevant and, locally, power 
laws are replaced by essential singularities. 

We consider two types of random extended surface perturbations associated with 
a pure bulk system: 

For the first type of perturbation, the probabilities p l and q L are constant and the 
couplings take the values: 

f' = {l < model#1 > < L2 > 

The perturbation has a vanishing average and the couplings decay towards their 
constant bulk value J. 

For the second type, the decay is associated with the probabilities whereas the 
amplitudes remain constant. The couplings follow the binary distribution: 

J l = J(l + /,a 2 ), /, = {J 9l =7i7-« ( model #2) (1.3) 

The average coupling takes the same form as in the HvL model and one recovers 
asymptotically a pure bulk system with couplings equal to J. 

The transverse field is always assumed to be constant, h t = h = 1, thus the bulk 
is critical at J = 1. 

In the present work, we study the surface magnetic properties in the marginal 
situation which occurs at k = | for model #1 and k = 1 for model #2 as well as for 
relevant perturbations corresponding to slower decays. 

The same type of problem with an aperiodic distribution of the couplings following 
some substitution sequence has been considered recently [17]. 

The structure of the paper is the following. The relevance of random extended 
surface perturbations is discussed in section 2. In section 3, analytical results about 
the surface magnetization of the inhomogeneous TIM are presented. The two models 
of inhomogeneous surface disorder are studied in section 4 for marginal perturbations 
and section 5 for relevant ones. The results are discussed in section 6. 

2. Relevance-irrelevance criterion 

We begin with a discussion of the relevance of random extended surface 
perturbations in a <i-dimensional classical system. The unperturbed system has a 



J, = J 



l + a 1 



(-1)" 
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free surface of dimension d — 1 at I = 0, perpendicular to the unit vector tt, and the 
position vector is written as r = lu + ry. The layered perturbation, which couples to 
the operator density 0(r) with bulk scaling dimension x^, is written as: 

- = J>(0*(r) , (2.1) 

r 

where (3 = (/c B T) _1 . The perturbation amplitude A(l) in the different layers are 
independent random variables such that: 

MOU = , [A(l)A(l')] av = ^ (model #1) 

A A- A- (2 - 2) 

[A(l)] av = , [A(l)A(lX v = S w + (1 - *„,) (model #2) 

Thermodynamic perturbation theory, up to second order, gives the following correction 
to the free energy: 

-pAF = -(3(V) + \fl 2 [(V 2 ) - (V) 2 ] 



= 5>(O(0(r)> + W A V)Ml')G H (ry) 



(2-3) 

( \ \U\AW\G ir.r') . 
r r,r' 

where (• • •) denotes a thermal average and is the connected two-point correlation 
function of the operator 0(r), both for the unperturbed system. 

In the case of model #1, the average of the first-order correction vanishes. The 
relevance of the perturbation is linked to the scaling behaviour of the amplitude in 
the average of the second-order correction: 

-piAFiX^^-EEw 8 ^' <2 ' 4) 

I l-ll 

where L d ~ l is the surface of the layers. The average perturbation density has 

dimension 2d — 1 wheras the correlation function has dimension 2x^ thus, under 
a length rescaling by a factor b = L/L' , 

A' 2 A 2 A 2 

A 1 _ h 2d-l-2 X4 , A l _ h 2y+-\ A l (n c\ 

and the second order perturbation amplitude transforms as: 

A'l = b 2y *- 2K - x A\, (2.6) 

where = d — x^ is the bulk scaling dimension of the variable conjugate to cf>. It 
follows that the perturbation is marginal when k = k x = y^ — |. It becomes relevant 
(irrelevant) when k < (>)k 1 - 

In the case of model #2, the average of the first-order correction is given by: 

-/?[AF«] av = £^(0(r)>. (2.7) 

r 
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The average perturbation density has now dimension d so that 

(2.8) 



(2.9) 

The perturbation is marginal when k = k 2 = y^, relevant (irrelevant) when k < (>)k 2 - 
One may notice that the average perturbation is of the HvL form, thus the relevance 
criterion is the same for both models to first order in the perturbation amplitude. A 
scaling analysis of the second-order correction contributed by the off-diagonal term in 
[A(l)A(l')] av leads to the same conclusion. Comparing the form of [A(l)A(l')] av for 
both models, one concludes that the diagonal second-order correction in model #2 
scales as A\ in equation (2.6) with 2k replaced by 1.6., clS 

[^diag = ^"^M^diag- (2-10) 

It follows that the second-order diagonal correction is relevant when k < 2y^ — 1. It 
is more dangerous than the first-order correction when > 1. 

In the case of the 2d Ising model with a thermal perturbation, such that (p(r) is 
the energy density cr r cr r+u , one has y^ = 1/v = 1 and marginal behaviour is obtained 
for k x = \ and k 2 = 1 as indicated in the introduction. 

The Id TIM, defined in equation (1.1) together with (1.2) or (1.3), which 
corresponds to the extreme anisotropic limit of the 2d classical systems discussed 
here [18], belongs to the same universality class, i.e., the relevance of the perturbation 
is the same for both problems. 



A' A 

2 = b d ~ x 4-fll 



and: 



A' 2 = b y *~ K A 2 



3. Surface magnetization and surface correlations 

Let us consider the imaginary time spin-spin autocorrelation function of the TIM: 

0,(r) = (0|of (r)of (0)|0) = K^f I0)| 2 exp[-r(^ - E )} , (3.1) 

i 

where |0) and \i) denotes the ground state and the ith excited state of 7i in equation 
(1.1), with energies E Q and respectively. With symmetry breaking boundary 
conditions, i.e., fixing the spin at the other end of the chain, such that a z L = ±1, 
the ground state of the system is degenerate and in the large-r limit we have 
lim,^^ Q^t) = m 2 , with a local magnetization given by: 

mi = (l\af\0). (3.2) 

To calculate m l we use the standard methods of Lieb, Schultz and Mattis [19] and 
transform 7i into a free-fermion Hamiltonian, 

^ = E e ,(^-^) > ( 3 - 3 ) 
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in terms of the fermion creation (annihilation) operators n q (i] q ). On a chain with 
length L and free boundary conditions, the non-negative fermion excitations e q satisfy 
the set of equations 



(3.4) 



= -h^ q {l) - Ji<P q (l + 1) , 

W) = - J i-i%( 1 - 1) - W) , 

with J = J L = 0. Introducing the 2L-dimensional vector V q with components 

U g (2/-l) = -^(/), ^(20 = ^(0, (3-5) 
the relations in equation (3.4) lead to an eigenvalue problem for the tridiagonal matrix: 



h x J x 
J, 



V 



\ 



J 



L-1 



h 



h, 



L 

/ 



(3.6) 



Taking the square of T, odd and even components of V decouple and one recovers 
two separate eigenvalue problems for ip and if). According to equation (3.4), changing 
<p q into — <p q in V , the eigenvector corresponding to — e q is obtained. Thus all the 
needed information is contained in that part of the spectrum of T with e q > 0. Later 
on we shall restrict ourselves to this sector. 

Fixing the surface spin at I = L amounts to put h L = in T, thus the lowest 
excitation e 1 = and the ground state of the system is degenerate. 

Using Wick's theorem, one may show that the local magnetization m z , with a 
fixed spin at / = L, is given by a determinant [8]. For the surface spin, at / = 1, the 
expression simplifies considerably and the surface magnetization takes the form [20]: 



m s (L) = m x = ^(1) 



L—l I /u 

i+ sn(^ 

i =lk=1 \ J k 



(3.7) 



This can be rewritten under the equivalent form 



L-1 



i=i 



(3.8) 



where fnf is the surface magnetization at / = L for the dual chain (for which the fields 
and the couplings are exchanged, h l <-> J t ) with a fixed spin now at / = 1 

We shall illustrate the usefulness of the last expression on the example of the HvL 
model which corresponds to the Hamiltonian (1.1) with h t = 1 and J L = J(l + al~ K ). 
In the marginal situation, k = 1, the product in (3.8) asymptotically behaves as 
YliLi ~ L a at the critical point, J c = 1. Furthermore, when L ^> 1, the couplings 
at the other end of the chain are asymptotically unperturbed, thus one expects that 
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mf(L) ~ L 2 , as for the homogeneous model at the ordinary transition. Then 
from (3.8) we deduce the continuously varying dimension 

which governs the scaling behaviour of the surface magnetization, m s (L) ~ L _Xm . 
This expression is valid only for a < \ since must be non-negative. When 
a > a c = x^ = 0, i.e., the critical surface magnetization remains finite for the 
semi-infinite system. It follows that our assumption about the scaling behaviour of 
should be corrected in the regime of first-order transition (a > \). It has to compensate 
the contribution of the product in (3.8), thus ra^ ~ L~ a . Since we are considering 
the magnetization of the dual chain, with couplings Jf(a) = l/J L (a) ~ J ; (— a), this 
means that when the couplings are sufficiently weak (a < — |) near / = 1, the finite- 
size-scaling behaviour of the critical magnetization at I = L is anomalous. 

When the perturbation is relevant, k < 1, the product of the couplings leads to a 
stretched exponential behaviour, 

m.-exp^^L 1 -^ , (3.10) 

which is clearly valid only when a < 0. Otherwise the dual magnetization term 
contributes, giving a size-independent result which corresponds to an ordered surface 
when a > [20]. 



4. Marginal extended surface disorder 

In this section we consider marginal random extended perturbations at the surface 
of the TIM. The transverse field is uniform, h l = h = 1, and the distributions of the 
couplings J l are given by either (1.2) with k = \ or (1.3) with k, = 1 for marginal 
behaviour. 

The finite-size-scaling behaviour of the average critical surface magnetization at 
J = 1 is deduced for both models from equation (3.8) for L 3> 1. As for the HvL model, 
one expects on the unperturbed end to scale as L~2 ? as long as the transition is 
continuous, while the average of the product behaves as [YliLi ^] av = Ilz^liWav 

For model #1, [J z ] av = J = 1 and [m s ] av ~ L"3, whereas nf=i[^] a v ~ L " 2 
for model #2. Thus we obtain the following values for the scaling dimension of the 
average surface magnetization: 

Klav=^ (model #1), { l i^]Z = 2 ' a * Ztl/2 (m ° del #2) • (ZL1) 

As for the HvL model there is a finite average surface magnetization at the critical 
point for strong enough local couplings, i.e., for a 2 > \ where lini L ^ 00 [m s ] av > 0. In 
this first-order regime the dual magnetization no longer scales as L -1 / 2 . Actually the 
correlation between the two factors on the r.h.s. of equation (3.8) modifies the scaling 
behaviour of the average of the product, leading to [i^] av = 0. One may notice that 
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Figure 1. Scaling dimension of the average surface magnetization as a function 
of the perturbation amplitude for model #1. The dotted line corresponds to the 
analytical result in equation (4.1). 




a 2 

Figure 2. Scaling dimension of the average surface magnetization as a function 
of the perturbation amplitude for model #2. The dotted line corresponds to the 
analytical result in equation (4.1). 



since the bulk Ising system has a correlation length exponent v = 1, equation (4.1) 
also gives the surface magnetic exponent [/? s ] av = [i s m ] ar 

The results for the average surface magnetization have been checked numerically, 
considering 10 6 realizations of randomly inhomogeneous quantum spin chains and 
computing the critical magnetization m s (L) given by (3.7) for chain sizes L = 2 5 to 2 14 . 
The exponents were deduced from an extrapolation of two-point approximants using 
the BST algorithm [21]. The variations of the average surface magnetic exponents, 
shown in figures 1 and 2, are in excellent agreement with the expected ones. 

Next we are going to study the probability distribution of the surface magnetization 
for the two models. Again we start from equation (3.8) in which, for large L, fn~, which 
refers to the unperturbed end, is expected to scale as L~z as long as the transition on 
the other surface is second order. Thus, in the second-order regime, the central limit 
theorem can be applied to the random variable lnm s + \ InL = Y^i=\ m (<^z) an d him s 
follows a normal distribution: 



P L (lnm s ) ~ 



a/2"7t var(lnm s ) 



exp 



(lnm s - [lnm s ] £ 
2 var(ln m s ) 



(4.2) 
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Figure 3. Probability distribution of lnm s calculated on 10 6 samples for model #1. 
The perturbation amplitude is a\ = 0.5 and the chain sizes are L = 2 6 , 2 8 , 2 10 , 2 12 
and 2 14 from right to left. The dotted line corresponds to the normal distribution in 
equation (4.2). 



where var(lnm s ) = [(lnm s — [lnm s ] av ) 2 ] av . This expression is valid for L ^> 1 and 
lnm s <C 0. It does not correctly describe the rare events, governing the average 
behaviour, which correspond to the tail of the distribution near lnm s = 0. 

The parameters of the distribution are easily deduced from equations (1.2) 
and (1.3) and read: 

1 L a 2 1 

[lnm s ] av ~ --lnL-^-i ~ --(l + a 2 )lnL, 



1=1 



(model #1) 



var(lnm s ) ~ a 2 lnL 



[lnm s 



ln(l + a 2 ) 



i=i 



I 



ln(l + a 2 ) 



InL. 



(4.3) 
(model #2) 



^ var(lnm s ) ~ ln 2 (l + a 2 ) InL 



The distributions P L (lnm g ), obtained through numerical calculations on 10 6 
samples for different chain sizes, are compared to the analytical expressions in figures 3 
and 4. The agreement is quite good at large sizes, although the normal distribution 
is strongly modulated for model #2. As shown in figure 5, the modulation involves 
a periodic function of period 1 in the variable lnm s /| ln(l + a 2 )|, where 1 + a 2 is the 
ratio of the two couplings in equation (1.3). The oscillations are log-periodic in the 
variable m s (see reference [22] for a recent review about log-periodic phenomena). 



The probability distribution is centered on [lnm s ] av = ln[m 



sJtyp 



InL, 



where [w s ] typ is the most probable or typical value of the surface magnetization. 
Consequently, the typical magnetization scales as a power of L, [ 

^"sJtyp ~ yP ' 

with a continuously varying typical exponent given by: 

)(l + a 2 1 ) (model #1), 



\x 



mJtyp 



X 



X 



mJ typ 
s 1 

mJtyp 



l/2-ln(l + a 2 ) a 2 <v^-l 
a 2 > y/e — 1 



(4.4) 



(model #2) 
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Figure 4. Probability distribution of In m s calculated on 10 6 samples for model 
#2 in the second-order regime. The perturbation amplitude is a,2 = —0.5 and the 
chain sizes are L = 2 6 and 2 14 from right to left. The dotted line corresponds to 
the normal distribution in equation (4.2). The actual distribution displays strong 
periodic oscillations around the normal distribution. 
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Figure 5. Variation of the period of oscillations in the distribution Pi(lnrn s ) for 
model #2 as a function of the marginal amplitude a^- The points give the periods 
deduced from the numerical data for the largest size, L = 2 14 , and the dotted line 
corresponds to | ln(l + a,2)\- 



In model #1, the typical surface magnetization vanishes continuously at the bulk 
critical point for any value of the perturbation amplitude. In model #2, for strong 
enough local couplings, a 2 > [o^ckyp = — 1 — 0.6487, there is a first-order surface 
transition, the typical surface magnetization remaining nonvanishing at criticality. In 
the range \ < a 2 < y/e — 1, at the bulk critical point, [m s ] av > whereas [m s ] typ = 
when L — > oo, a behaviour which is consistent with the fact that [lnm s ] av < ln[m s ] av . 

The numerical results for the scaling dimension of the typical surface magnetization 
are shown in figures 6 and 7. Two-point finite-size approximants were extrapolated 
using the BST algorithm and we used the same samples as for the average 
magnetization. The agreement with the expressions given in (4.4) is quite good. 
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Figure 6. Scaling dimension of the typical surface magnetization as a function of 
the perturbation amplitude a± for model #1. The dotted line corresponds to the 
analytical result in equation (4.4). 
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Figure 7. Scaling dimension of the typical surface magnetization as a function of the 
perturbation amplitude for model #2. The dotted line corresponds to the analytical 
result in equation (4.4). 



5. Relevant extended surface disorder 

Let us now examine the modified surface critical behaviour induced by relevant 
extended surface disorder which, according to the disussion of section 2, corresponds 
to a decay exponent k in equations (1.2) and (1.3) lower than k x = 1/2 for model #1 
and k 2 = 1 for model #2, respectively. 

Like in the marginal case, the average surface magnetization obtained by averaging 
equation (3.8) would scale with an exponent [a^Jav = 1/2 fc> r model #1 if the dual 
magnetization was unaffected by the perturbation. A numerical finite-size scaling 
study with a 1 = 0.5 and k = 0.25, shows that the actual behaviour is a stretched 
exponential one since ln[m s ] av ~ L 1 / 2 , as shown in figure 8. Thus the slowly decaying 
perturbation changes the critical behaviour on the second surface. By symmetry, the 
behaviour is the same for negative values of a 1 . Numerical studies for other values of 
k indicate that the argument of the stretched exponential involves the power L 1_2k . 

For model #2, taking only into account the product of the couplings, one obtains: 
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Figure 8. Finite-size scaling of the average (dot-dashed lines) and typical (dashed 
lines) surface magnetization with relevant perturbations illustrating the stretched 
exponential behaviour. The parameters are a\ = 0.5 and k = 0.25 for model #1 
(circles), ai = —0.5 and k = 0.5 for model #2 (squares). Averages were taken over 
10 6 samples. 
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Figure 9. Probability distribution of lnm s calculated on 10 5 samples for model #1 
with a relevant perturbation. The amplitude is ai = 0.5 and the decay exponent 
k = 0.25. The chain sizes are L = 2 6 , 2 10 and 2 14 from right to left. The dotted 
lines correspond to the normal distribution in equation (4.2) when the contribution 
of the dual magnetization to In m s is neglected. 



exp 



l \-k 



K 



(5.1) 



Hence the average surface magnetization has a stretched exponential dependence on 
L for negative values of the amplitude a 2 . The numerical data are in agreement with 
this behaviour as shown in figure 8 for a 2 = —0.5 and k = 0.5. The value of the slope, 
—0.99 ± 0.04, is close to the analytical result, —1. Thus the dual magnetization is 
likely to contribute here only through a power of L in front of the exponential. 

When a 2 > 0, the divergence of [w s ] av with the system size signals that the dual 
magnetization term in (3.8) can no longer be neglected. The surface transition is then 
first order. This was checked numerically at a 2 = 0.5 where the average critical surface 
magnetization is indeed size-independent at large L values. Contrary to the case of a 
marginal perturbation where a mininum enhancement of the couplings was needed to 
maintain surface order at the bulk critical point, here a first-order surface transition 
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Figure 10. Probability distribution of In m s calculated on 10 5 samples for model #2 
with a relevant perturbation. The amplitude is <Z2 = —0.5 and the decay exponent 
K = 0.5. The chain sizes are L = 2 10 , 2 12 and 2 14 from right to left. The dotted 
lines correspond to the normal distribution in equation (4.2) when the contribution 
of the dual magnetization to lnm s is neglected. 



occurs for any positive value of the perturbation amplitude. The same behaviour is 
obtained with the HvL model [20] . 

The probability distribution P L (lnm s ) can be obtained in the same way as in the 
marginal situation when the transition is continuous, provided that lnm s is dominated 
by the contribution of the product in (3.8). Then one recovers the normal distribution 
of equation (4.2) with the following parameters: 



^2 2 



1-2k 



var(lnm s ) ~ 



a i 7-1-2K 



s/ 1 2k 

L 



1 = 1 



n n In (1 + a 2 ) r i_ K 
var(lnmj ~ ; — L K 

S/ 1-K 



(model #1) 

(5.2) 

(model #2) 



The numerical data, shown in figures 9 and 10, were deduced from 10 5 realizations 
and, as above, with a 1 = 0.5, k = 0.25 for model #1 and a 2 = —0.5, k = 0.5 for model 
#2. 

In the case of model #1, the shift with respect to the analytical results confirms a 
weak relevant contribution of the dual magnetization term in (3.8). This is confirmed 
in figure 8 where the slope of ln[m s ] t versus L 1_2k is 0.266 ± 0.003 to be compared 
to the value —0.25 given by (5.2). 

The agreement is much better for model #2. Here the slope of ln[m s ] t versus 
L l ~ K in figure 8, —1.401 ± 0.002, has a smaller relative deviation from the expected 
value, —1.386. One may notice that the probability distribution no longer displays 
the strong log-periodic modulations observed in the marginal case. 
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6. Discussion 



The two models of extended surface disorder studied in this paper lead to a 
surface critical behaviour which is quite similar to that obtained in the HvL model, 
with continuously varying exponents in the marginal case and a stretched exponential 
behaviour for relevant perturbations. The main differences are the non-self-averaging 
behaviour of the random models and the the log-periodicy obtained with model #2. 

This log-periodic behaviour is easily understood if one considers the expression of 
the surface magnetization in equation (3.8). At criticality h t = J = 1 and, according 
to (1.3), In J z is either for unperturbed couplings or ln(l + a 2 ) for perturbed ones. 
A sample a with n a modified couplings contributes to the probability distribution 
P L (lnm s ) at a value n a ln(l + a 2 ) + lnm^ a of lnm s . If the fluctuations of the last 
term from sample to sample are sufficiently weak compared to the interval ln(l + a 2 ) 
a modulation with period ln(l + a 2 ) is obtained for P L (lnm s ). The argument does 
not hold for model #1 since the values of the couplings vary from site to site. The 
absence of modulation for model #2 with a relevant perturbation gives a lower bound 
for the influence of the disorder on the second surface. 

The connection to the HvL model can be clarified through the introduction of 
effective HvL models describing either the average or the typical behaviour. As long 
as the size dependence of the dual magnetization in equation (3.8) is the same as for 
an unperturbed surface, one may introduce effective HvL interactions such that: 

Wav = Wav , Myp = «P(MD • (6.1) 

Equations (1.2) and (1.3) then lead to: 



[ J l\typ — J (l 2 /2k) 



WIS = J (i + £) . 



[«^]typ — ^ 



1 + 



/' 

hi(l + a 2 ) 



(model #1) 

(6.2) 

(model#2) 



For model #1, k HvL = 2k and the random perturbation is marginal at k 1 = 1/2. 
The effective amplitude vanishes for averaged quantities whereas it is always negative, 
with [a]ty P = — «i/2, for typical quantities. Hence the surface transition is always 
continuous and the marginal exponents in (4.1) and (4.4) follow from (3.9) with the 
appropriate effective amplitude. 

For model #2, k HvL = k leads to k 2 = 1. The effective amplitudes are [a]|f = a 2 
and Htyp = m (l + a 2)- The continuously varying exponent can be deduced from 
the HvL results when the transition is continuous, in this case even for relevant 
perturbations (compare (5.1), (5.2) and (3.10)). 

One may go further and conjecture the form of the scaling dimension of the surface 
energy density in the second-order regime. Its value x s e = 2(1 — a) for the HvL model 
[23] translates into: 
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= 2(1 -a 2 ) 
, = 2[l-ln(l + a 2 )] 



(model #1) , 



a 2 < 1/2 
a 2 < \/e — 1 



(model #2) . 



(6.3) 



As already mentioned, the inhomogeneously disordered Ising quantum chain 
corresponds to the extreme anisotropic limit of a 2d classical Ising model with 
correlated disorder along the rows parallel to its surface. It is easy to apply the 
considerations of section 2 to the case where the perturbations of the different bonds 
(or sites) at a distance / from the surface of a (i-dimensional classical system are 
independent random variables distributed according to (1.2) or (1.3), i.e., with d- 
dimensional inhomogeneous disorder. 

For model #2, to first order, the perturbation amplitude A 2 transforms like in (2.9) 
with a scaling dimension — k when it couples to the field 0(r) in d dimensions. The 
off-diagonal second-order correction has the same behaviour. The diagonal second- 
order correction, with scaling dimension 2y^ — n — d, cannot be more relevant than the 
first-order correction since y^ < d. Thus the relevance of the perturbation remains 
the same as for the HvL model. 

For model #1, however, the scaling dimension of the second-order perturbation 
amplitude A\ in equation (2.6) changes to y A 2 = 2y^ — 2k — d. Thus a thermal 
perturbation is marginal when k = \ jv — d/2 = a/(2v) where a is the bulk specific 
heat exponent. For the 2d Ising model, with a = 0, the perturbation is irrelevant 
as soon as k > and marginal for homogeneous disorder, a well-known result [24]. 
The 2d g-state Potts model would be more interesting since marginal behaviour then 
corresponds to k = 1/5 for q = 3 (is = 5/6) and k = 1/2 for q = 4 {y = 2/3) [25]. One 
could then interpolate between a marginal inhomogeneous perturbation and a relevant 
homogeneous one. Another case of interest is the 2d g-state Potts model with q > 4. 
Here the pure system has a first-order bulk transition [26,27], the surface transition 
is continuous [28, 29] and homogeneous disorder changes the bulk transition from first 
to second order [30-32]. With v = 1/2 at the discontinuity fixed point [33,34], the 
critical decay exponent for marginal behaviour is k = 1. Varying k from 1 to could 
lead to interesting new phenomena. 

To conclude, let us mention that it would be useful to adapt Fisher's 
renormalization group approach to the type of perturbations considered in this work. 
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